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Abstract 

The acceleration of charged particles in magnetized plasmas is considered during turbulent multi- 
island magnetic reconnection. The particle acceleration model is constructed for an ensemble of 
islands which produce adiabatic compression of the particles. The model takes into account the 
statistical fluctuations in the compression rate experienced by the particles during their transport 
in the acceleration region. The evolution of the particle distribution function is described as a 
simultaneous first and second-order Fermi acceleration process. While the efficiency of the first- 
order process is controlled by the average rate of compression, the second order process involves the 
variance in the compression rate. Moreover, the acceleration efficiency associated with the second- 
order process involves both the Eulerian properties of the compression field and the Lagrangian 
properties of the particles. The stochastic contribution to the acceleration is non-resonant and can 
dominate the systematic part in the case of a large variance in the compression rate. The model 
addresses the role of the second-order process, how the latter can be related to the large-scale 
turbulent transport of particles and explains some features of the numerical simulations of particle 
acceleration by multi-island contraction during magnetic reconnection. 

PACS numbers: 96.60.qe, 52.35.Vd, 52.65. Cc, 96.60.Iv 
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I. INTRODUCTION 



The production of non-thermal particles in magnetized plasmas is an ubiquitous complex 
phenomenon which is believed to involve also magnetic reconnection. Magnetic reconnection 
is the process that controls the conversion of magnetic energy into kinetic energy[l|; it is the 
driver of impulsive phenomena such as solar flares, substorms in the Earth's magnetosphere 
and disruptions in laboratory fusion devices. The relation between magnetic reconnection 
and particle acceleration has been extensively discussed in the terrestrial magnetosphere 
based on in-situ observations Moreover, X-ray observations and studies of the energy 

budget during solar flares indicate that a significant fraction of the magnetic energy released 
in a flare is carried by the accelerated 10 — 100 keV non-thermal electrons [s}]. However, how 
particles can be accelerated in large numbers to high energies as the magnetic field lines 
reconnect remain an outstanding problem. 

Numerical PIC simulations aiming to address the problem of particle acceleration dur- 
ing magnetic reconnection in a self-consistent manner, have confirmed that particles are 
efficiently accelerated in the vicinity of the X-line by reconnection electric fields . How- 
ever, an important limitation of X-type acceleration mechanisms is that they hardly explain 
alone the large number of accelerated particles, in particular during a solar flare, because 
the volume occupied by a current sheet where the strong electric field capable of particle 
acceleration is present, is quite small. 

Hence, intensive efforts have been made to understand the role played by the region 
inside the separatrices for particle acceleration, leaning toward the idea of 0-type accelera- 
tion mechanisms which take advantage of the closed geometry of the field within magnetic 
islands. Drake et al. [8^ have developed a model of particle acceleration which is based on 
the dynamical motion of the islands. They show that particles trapped in the contracting 
magnetic field of the islands are adiabatically compressed and therefore can be efficiently 
accelerated through a first-order Fermi process. In addition, many studies have revealed 

nn nn 

the importance of magnetohydrodynamics turbulence [9l-ll ll| and plasmoid dynamics |12l . Il3l | 
with regards to particle acceleration in a reconnecting plasma. Current sheets are naturally 
prone to tearing and their fragmentation lead to the formation of magnetic islands having a 
complex multiscale and intermittent dynamical behavior fl^ . Iisl . Our goal in this Letter is 
to study the effect on particle acceleration of the ensemble of contracting islands and develop 
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a simple model of particle acceleration during turbulent magnetic reconnection. 

When a magnetic island changes its length L at the velocity V = dL/dt particles that 
are trapped within the island change their speed v according to dv/dt = —{aV/L)v = 
Wv. This relation, derived in is a consequence of the conservation of the longitudinal 
action for particles trapped within an island and the coefficient a represents the relative 
magnitude of the reconnecting magnetic field, i.e. a = {SB/BoY. As a result, magnetic 
islands that are contracting at a speed of the order of the Alfven speed V ~ Va, accelerate 
the trapped particles through adiabatic compression, provided f ^ Va- The acceleration 
rate associated with this first-order Fermi process is given by aVa/Lo, where Lq is the 
typical length of the islands. If on the contrary the islands are expanding, then first-order 
adiabatic deceleration of the particles will result at the same energy independent rate. When 
contraction is magnetically favorable and when the energy gained from the magnetic field by 
the particles is balanced with energy losses, including transport losses and/or back-reaction 
of the accelerated particles, power-law distribution in particles energy may be obtained, 
which are determined by standard techniques, as was originally done in jsj. 

In the case of a first-order Fermi process, the rate of energy gained by the particles is 
proportional to the mean compression (W), where the brackets () denote an average over the 
ensemble of islands in the system, possibly weighted by the relative number of islands that are 
undergoing contraction mere existence of this average, or the range of possible 

contraction rates, suggests to consider also the effect of the finite variance in the adiabatic 
compression experienced by the particles in the sea of islands. Indeed, the contraction rate 
changes in time due to firehose condition 16|, so an assemble of contracting island will have 



non-zero variance. Further, PIC simulations 



12| emphasize the bouncing motion of merged 



islands, so that a contracting motion of an island is followed by an expanding motion. 

For an ensemble of multiple contracting islands, the presence of non-zero average (W) 
and non-zero {{W — (VT))^) leads to both first and second order accelerations. The mean 
controls the first order Fermi acceleration and additional statistical acceleration occurs at 
a rate proportional to the variance of the compression, also when the mean compression 
rate (W) is non-zero. A continuity equation can be written for the omnidirectional particle 
distribution function F{p,t) = Anp^ f{j),t) 



dt dp Wdt 



0, (1) 
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where p is the particle momentum with the time rate of change in momentum given by 

dp aV 



dt L 



-p = Wp. (2) 



In 



16| . a term modehng the effect of escape of particles out of the acceleration region is 
also included in Eq.(Il]). 

Let us consider the case where the compression rate is small and assume first that it is a 
function of time only with zero average, i.e. bouncing motion of the islands 

{W{t)) = 0, (3) 

and its correlation function decays exponentially, 

C{t) = {W\t))eM-t/rc). (4) 

In this cases, even when on average the islands are neither contracting nor expanding, 
i.e. iyV) = 0, there is a stochastic acceleration effect that remains operative. Although, 
the particles do not experience any systematic change in their energy, the average particle 
energy could still grow with the acceleration efficiency associated with a second-order Fermi 
process is proportional to the variance of the compression {W"^). Therefore, we obtain that 
the mean omnidirectional distribution function Fo{p,t) obeys the diffusion equation 

with the diffusion coefficient in momentum space given by 

POO 

D= / dtCit)=T,{W\t)). (6) 
Jo 

The mean distribution function Fo{p,t) solution of Equation is the normal distribution 
with respect to the variable u = ln(p/po)- Indeed, the particle dynamics is described by the 
Langevin equation du/dt = W(t) with {W(t)) = 0. Therefore, Fo{u,t) satisfies the standard 
diffusion equation 

dt du^ ' 

which also confirms that fluctuations in the compression rate are responsible for the growth 
of the variance in the momemtum distribution function. 

An account for the effect on the stochastic acceleration of the spatial transport of particles 
in the pulsation field of the islands may be given on the following basis. Let us shrink the 



volume of each island into a point, this point being characterized by its compression W{'x, t) 
with X being the position of the center of the islands. Moreover, we envisage a situation where 
the large scale spatial transport of particles in the volume filled by the islands is turbulent and 
diffusive. Therefore, the particle dynamics is modeled by the following Langevin equations: 

I^CW; (8) 

with (C(t)) = and (OWO(O) = 25ijKT 5{t - t'), {W{^,t)) = and {W{0,0)W{x,t)) = 
C(x, t), kt is the spatial diffusion coefficient. Here C(x, t) is the Eulerian correlation function 
associated with the compression/expansion field W{yi,t) of the islands which is supposed 
to be homogeneous and stationnary. The Eulerian correlation function depends on three 
parameters that characterize the statistics of the (isotropic) compression/expansion field : 
the variance (Vr^(x, t)) = C(0,0), the correlation time Tc, which is the decay time of the 
Eulerian correlation and the correlation length Ac, which is the decay length. So the particles 
have probability to stay within the island or escape. As noted [8], the gyration radius of the 
particle increases near the separatrix, which in turn increases the probability of a particle to 
escape the island. The Langevin equations ([H]) are doubly stochastic in the sense that both 
the position x(t) of the particles and the compression field W^(x, t) are stochastic processes. 

With the spatio-temporal statistics of the compression being specified via C(x, t), the 
problem is to calculate the diffusion coefficient in momentum space (when the latter exists) 
and to determine the form of the distribution function. The diffusion coefficient D is related 
to the time integral of the Lagrangian correlation function[l3], viz. 



D= / dtCL{t), (9) 
Jo 

where the Lagrangian correlation function Ci{t) is defined via 

C^(t) = (l^(0,0)l^(x(t),t)), (10) 

where x(t) is a solution of Eqs.([8]). The exact result ([9]) is a simple consequence of 
the definition D = {l/2)d{u^) / dt combined with the second equation in ([8]). Indeed, 
(^2) = jldt' f^dt"{W(t')W{t")) = 2j^dt'CL{t'){t-t') and letting t oo (when the 
integral converges) gives Eq.(l9]). It also follows from Eq.(j9]) that the diffusion coefficient in 
momentum-space can be expressed as 

D = rL{W'{^,t)), (11) 



where tl is the Lagrangian correlation time, i.e. the correlation time of the compres- 
sion/expansion field which is experienced by the particles along their trajectory. The prob- 
lem remains to connect Lagrangian and Eulerian statistics, i.e. to determine the functional 
dependence of with and Ac. To this purpose, let us write the Lagrangian correlation 
function f llUp in the equivalent form 

CLit) = y"c/x(Vr(0,0) Vr(x,t)5[x-x(t)]). (12) 

A relation between the Lagrangian correlation Crit) and the Eulerian correlation C(x, t) is 
obtained by invoking a procedure due to Corrsin 18, lol in which x(t) is replaced by its 
statistical average, so that we may replace 5[x — x(i(:)] in equation ( TT^ by ((5[x — x(t)]). This 
leads to the factorization Ci(t) = / d:x.{W{0,0)W{'x,t)) (5[x — x(t)]). Hence, an expression 
for the diffusion coefficient in momentum space is found which is given by 

D = J dt J dxC{x,t) P{x,t) , (13) 

where P(x, t) = (5[x — x(t)]) is the conditional probability for a particle to be under the 
influence of a magnetic island located at the position x at time t provided that this particle 
was at X = at t = 0. Equation (|T3l) shows that D is the integral of the product of 
two quantities: C(x, t), the Eulerian correlation function, which characterizes the statistical 
properties of the compression fleld W^(x, t), and the probability function P(x, t), describing 
the spatial transport of particles in the acceleration region. Here, P(x, t) is the solution of 
a standard diffusion equation with diffusion constant /ty, i.e. 

P(x,t) = ^^j^tr/2 ^M-\^\V^^Tt), (14) 

but the procedure can be generalized to more complex transport models. 

Let us further take the illustrative example of an isotropic correlation function of the 
form given by 

C(x,t) = (iy2(x,t))exp(-|x|VA^-t/r,). (15) 
From Eq. f|T3l) . we obtain that 

D = {W\^,t)) r dtexp{-t/r,)il + ^^)-'/\ (16) 

Jo 

Therefore, in the weak spatial diffusion limit, where kt ^ A^/4rc, the momemtum diffusion 
coefficient is given by 

D^T,{W\^,t)). (17) 
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FIG. 1: Particle distribution function (top panel) and the spectral index of the distribution (bottom 
panel). The solutions of Eq.© for Dt = 3 (solid line), Dt = 1 (dashed line), and Dt = 0.5 (dash- 
dotted line). All distributions are normalised so that f F{){p^t)dp = 1. 

This is the case already given by Eq.(l6]) corresponding to Eq. lfTTI) with ~ Tc- However, 
in the opposite, strong spatial diffusion limit, where k^t ^ -^c/^'^c, then 

^~^(W^'(x,t)), (18) 

corresponding to the Lagrangian correlation time being of the order of the spatial transport 
time scale, i.e. ~ A^/k^. In this strong spatial diffusion limit, the stochastic acceleration 
efficiency is governed both by the Eulerian properties of the compression field and the 
Lagrangian properties of the particles. 

Let us notice that the diffusion coefficient in momentum space may also be expressed 
as D = // dkduj S{\i.,uj) nrk"^ /[^"^ + (^t^^)^], where S'(k,co') is the spectrum of Vr(x,t), 
i.e. the Fourier transform of the correlation function C(x, t). It can be clearly seen from 
this expression for D that the integral may diverge for scale-free power-law spectra such 
as S{k,u) cx k~'^6{u). This is the signal that the turbulent acceleration process cannot be 




described as a standard diffusion in m— space as in Eq.Q. This situation has been dubbed 
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2l| . Here, we focus on the second-order process 



Fermi acceleration of fractional order in 
with D finite. 

The statistical effect discussed above can be felt also in addition to the systematic energy 
change. Indeed, when both the mean and the variance of the compression are finite, the 
first and second-order Fermi processes operate together. In this case, Fo{u,t) obeys an 
advection-diffusion equation in velocity space, 

dFo(u,t) d ^ , . d'^ ^ , . . 
^^ + «i^^oKt) = a2^Fo(n,t), (19) 

where u = ln(p/j9o) and where the coefficients of systematic and stochastic acceleration are 
given by 

ai = (iy(x,t));a2 = ri((lV(x,t) - {W{^,t))f), (20) 

respectively. When the islands contract on average, the distribution function Fo{u,t) shifts 

toward large m at a rate given by ai while the variance of FQ{u,t) grows at a rate given by 

02 and the stochastic component to the acceleration process dominates the systematic part 

when a2 ^ fli. The time-dependent solution Fo(u,t) of the advection-diffusion equation 

( IT9|) is the normal distribution in the variable u — ait. 

Although the time-dependent solution is not a power law, but only asymptotically at 

t — 7- oo, the characteristic solutions and spectral indices 

_ d\nFo{p,t) _ Inp/po 

dlnp Dt ' ^^^^ 

are found for a few values of Dt and are presented in Figu re [ll The values appear to be 



similar to those obtained in numerical simulations e.g 
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values in solar flares [5j than for example in 11|. 



16| and closer to the observed 



In summary, we show that both the first and second-order Fermi acceleration process can 
operate together to increase the particle energy when the acceleration region consists of an 
ensemble of contracting islands. In the case when islands are both contracting and expanding 
with zero mean effect, only the second-order acceleration process operates. However, even 
when contraction is dominant, the second order effect can be substantial. The stochastic 
component to the acceleration corresponds to a non-resonant mechanism according to the 
classification scheme established in [21]. It involves the turbulent transport properties of the 
particles in the acceleration region and becomes more efficient for higher levels of variance 
in the compression rate. 
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